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ABSTRACT 

We show that the minimally coupled massless scalar wave equation in the background of 
an six-dimensional extremal dyonic string (or D1-D5 brane intersection) is exactly solvable, 
in terms of Mathieu functions. Using this fact, we calculate absorption probabilities for 
these scalar waves, and present the explicit results for the first few low energy corrections 
to the leading-order expressions. For a specific tuning of the dyonic charges one can reach 
a domain where the low energy absorption probability goes to zero with inverse powers 
of the logarithm of the energy. This is a dividing domain between the regime where the 
low energy absorption probability approaches zero with positive powers of energy and the 
regime where the probability is an oscillatory function of the logarithm of the energy. By 
the conjectured AdS/CFT correspondence, these results shed novel light on the strongly 
coupled two-dimensional field theory away from its infrared conformally invariant fixed 
point (the strongly coupled "non-critical" string). 
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1 Introduction 



One of the more intriguing recent developments in string theory and M-theory has been the 
conjecture that when the theory is placed in a background that is of the form of a product 
of an m-dimensional anti-de Sitter AdS m spacetime and an re-sphere S n , it is dual to a 
superconformal theory on the boundary of the AdS m spacetime |IJ. The AdS m x S n sphere 
structure arises in the near-horizon limit of non-dilatonic extremal p-branes 0, |3| , and thus 
one arrives at a conjectured relation between aspects of the p-brane supergravity solution 
and the superconformal theory on the AdS m boundary. 

This conjectured relationship has been explored in a number of contexts. One partic- 
ular aspect has been probed by studying the absorption probabilities and scattering cross- 
sections for massless fields propagating in the supergravity background describing p-brane 
solitons J|, ||, In particular, one is principally interested in studying the absorption 
probability for incident waves with low frequency u. For example, the leading low-energy 
absorption probability for t = partial wave of the the minimally coupled massless scalar 
field (dilaton) determines its two-point correlation function at the conformally invariant 
infrared fixed point of the strongly coupled (m — l)-dimensional field theory 0, ||, [| 1C]. 
Sub-leading corrections are expected to shed light on the effects of (irrelevant) perturbations 
away from the conformally invariant fixed point. 

In many of the examples that have been studied, the wave equations turn out not 
to be explicitly solvable in closed form, and as a result various techniques for matching 
solutions in an overlap between approximately-solvable inner and outer regions have been 
used 0, 0, HI, 0, H, |, |, H, |0], 0, H, H, |J]. The possibility of making 

such overlapping approximations arises because of the low-frequency nature of the incident 
waves. One can then obtain expressions for absorption probabilities as power-series in oj. 
In most cases where such approximations must be made, it is difficult to obtain more than 
a few sub-leading corrections to the leading-order absorption probabilities. 

On the other hand, in certain special cases (of extremal »-brane systems) the wave 
equations turn out to be exactly solvable, and in these situations one has more complete 
control over the calculation of corrections to the leading-order absorption probability. One 
particularly interesting example arises for the extremal D3-brane, where the wave equation 



for massless scalars turns out to be precisely the (modified) Mathieu equation [25]. This 
equation has been much studied by mathematicians, and techniques have been developed 
that are directly adaptable to the problem of calculating the absorption probability. Indeed, 
the work of p(l formed the basis of the approach used in p5| for calculating the D3-brane 
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absorption probability. 

In this paper we consider another example where the wave equation is exactly solv- 



able, namely the extremal dyonic string [28, 29] in six dimensions. This configuration can 
alternatively be viewed, by diagonal dimensional reduction, as being equivalent to an ex- 
tremal 2-charge black hole in D = 5, a Dl-brane/D5-brane intersection in D = 10, or an 
M2-brane/M5-brane intersection in D = 11 J27j. From the point of view of the AdS/CFT 
conjecture this is an especially interesting case, since the near-horizon geometry of the dy- 
onic string approaches AdS3 x S 3 , and thus the superconformal field theory [^] is defined 
on the two-dimensional boundary of AdS3. A controlled method to calculate the absorp- 
tion probability to an arbitrary order in energy corrections is thus of interest; such a result 
would provide information on the strongly coupled two-dimensional field theory perturbed 
away from the infrared conformally-invariant fixed point by (irrelevant) operators, i.e. a 
strongly-coupled "non-critical string" would be probed. 

Low energy absorption probabilities for the Dl-brane/D5-brane system, or equivalently 
the six-dimenesional dyonic string, were considered in |2^] , by using the approximate meth- 
ods described above, which involve a division of the spacetime into overlapping inner and 
outer regions. In particular, the next-to-leading order energy correction to the s-wave ab- 
sorption probability was given there. 

In this paper we show that the wave equation for the minimally coupled massless scalar 
in the extremal six-dimensional dyonic string background has in fact an exact solution; a 
straightforward change of variables transforms it into the modified Mathieu equation. Thus 



we can again, as was done for the D3-brane [25|, make use of the results in |26|], and calculate 
the absorption probabilities to an arbitrary order in the energy expansion. However, there 
are some significant differences between the nature of the low energy expansion of the results 
presented in this paper and those of the extremal D3-brane, which we shall explain later in 
the paper. 

The paper is organised as follows. In section 2, we write down the scalar wave equation 
in the background of an extremal dyonic string, and show that it can be transformed into the 
Mathieu equation. In section 3, we discuss the details of how to solve the Mathieu equation, 
paying particular attention to some subtleties associated with the particular parameter 
values that arise in the dyonic background, and their energy dependence. We obtain the 
absorption probability as an energy expansion for a fixed dyonic string configuration. In 
section 4, we explore some different regions in the parameter space, for which the small 
expansion parameter is governed by the ratio of the two dyonic string charges. By tuning 
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the ratio of the two charges one probes regimes whose absorption probabilities exhibit novel 
functional dependence on the energy parameter. 



2 Scalar wave equation for the dyonic string 

We shall consider a minimally coupled scalar field propagating in the background of an 
extremal dyonic string in six dimensions, for which the metric is 

dsl = (H e Hm)" 1 / 2 (-dt 2 + dx 2 ) + (H e H m ) 1 / 2 (dr 2 + r 2 dQ 2 3 ) , (2.1) 

where 

H e = l + ^, Hm = l + ^. (2.2) 
The wave equation □<!> = can be separated, by assuming coordinate dependence 

$(t,r,0i) = <f>{r)Y{e i )e- i " t , (2.3) 

where Y(9j) denotes an harmonic on 5 3 , satisfying V 2 Y = —£(£ + 2) Y on the unit 3-sphere. 
The radial function <p(r) therefore satisfies the equation 

d?cj> 3 # / , l(l + 2) x 

-T)+--r+ w H e H m ^ — = 0. (2.4) 

ar z r dr \ r z t 

Defining p = ojr, the equation becomes 

g + »g + ( 1 + £±£4£g±2 + £ft),- , (2.5) 
dp L p dp ^ p z P ' 

where Ag = u> 2 Q e and = lo 2 Q m . Note that we would obtain exactly the same equation 
(|2.5| ) if we were to consider the intersection of a Dl-brane and a D5-brane in D = 10, or 
for such an intersection with NS-NS branes, or an intersection of an M2-brane and an M5- 
brane in D = 11. We would also encounter the same equation if we considered a 2-charge 
black hole in D = 5. These equivalences hold because the associated brane configurations 
are related to one another by diagonal dimensional reduction, and the wave equation is 
invariant under such reductions on the world volume. 



It is easy to see that the equation (2^5) can be transformed into the Mathieu equation, 
by defining 

<Kp) = -*(p), P = VAe X m e~ z , (2.6) 
P 

whence we obtain 

V + (8A 2 cosh(2z) - 4a 2 ) * = . (2.7) 
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Here, we have defined 



a 2 



X 2 = jX e X m = \lo 2 \J Q e Q m , (2-8) 



Ae A m / Q e / Qm Q e + Q r 



Am A e y Qm y Qe vQe Qm 

Note that A > 2, and that A will be treated as a small parameter (corresponding to low- 
energy waves). 

The result that we have found here is in many respects similar to the one for the D3- 



brane obtained in [25|, where it was shown that the scalar wave equation in the D3-brane 
background could be reduced to the Mathieu equation. However, a significant difference 
arises in our case, due to the fact that now the parameter a, appearing in the Mathieu 
equation, itself depends on the small (energy dependent) parameter A. Although we shall 
be able to follow the same general strategy for solving the equation as the one described in 
p6| , |25|| , some adjustments will be required in order to take into account the additional A 
(energy) dependence of parameter a. 

The solutions to the Mathieu equation are constructed as an iterative sequence of terms 
which correspond to a power series in the small parameter A. There are various ways in 
which one can view the relation between a and A in equation ( |2.8| ), The simplest of these 
is to think of it as defining a as a power series in A, so that one has a = \(l + 1) + 0{X 2 ). 
In this viewpoint, one ensures that A is small by choosing the frequency lo to be sufficiently 
small, with no constraints on the charges. We calculate the absorption probabilities within 
this framework in section 3. One can alternatively think of a as being a freely-specifiable 
parameter (independent of A), where one ensures that A, now viewed as being defined by the 
relation between a and A in equation (2.8), is small, by choosing the charges so that A is 



sufficiently large. We calculate the absorption probabilities within this framework in section 
4. (In order to distinguish the somewhat different interpretations of the A parameter in the 
two viewpoints, we replace A by A in section 4 when we are viewing it as a derived quantity.) 
It is of particular interest is to explore the limit where a — > 0, since this approaches the 
boundary of the region a 2 < 0, which has a novel low-energy behaviour for the absorption 
probability. 
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3 Solving the Mathieu equation 



A procedure for solving the Mathieu equation was developed in [26|, and formed the basis of 
the approach used in [p5[ . In order for us to present our notation and results, and to explain 
certain subtleties that arise, it will be necessary for us to review the procedures developed 
in [26, pq] . At the same time we shall elaborate on some delicacies in the calculations that 
were not explicitly discussed in those papers. 

3.1 The Floquet expansion 

The Mathieu equation can be solved by mapping the problem into one of solving a difference 
equation for coefficients in a series expansion. Specifically, the solution can be expressed in 
the Floquet form 

oo 

^(z)= Y, (-l) n C(n + / u)e 2(n+ ^ ) % (3.1) 

n=— oo 

where \x is a certain constant, related somewhat non-trivially to the parameters in the 
Mathieu equation, which will be determined later. The quantity 2/x is known as the Floquet 
exponent, or Mathieu characteristic exponent. Substituting ( |3.1| ) into ( |2.7D , one finds that 
the coefficients C{n + \x) satisfy the recursion relation 

C{x + 1) + C(x-1)= C(x) . (3.2) 



Convergence of the series (|3.1|) requires that C(x + 1) << C(x — 1), and hence an approx- 
imation C (0) (x) to the solution of this three-term recursion relation can be obtained by 
solving the two-term relation 



C<°> (x - 1) = ^ " 2 C (0) (x) . (3.3) 
A 2 



This is easily seen to lead to 



\2x 

C w (x) = — ^— . (3.4) 

W T(x + a + l)T(x-a + l) y J 



A solution to the full recursion relation (|3.2| ) can now be obtained by writing 

C(x) = C m (x)B(x) , (3.5) 



where, substituting into ( |3.2| ), one finds that the coefficients B{x) must satisfy 

B {x) - B(x + 1) = - ^B(x + 2). (3.6) 

{{x + iy — a z ) ({x + 2y — a z ) 
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Since A is treated as a small parameter, we can therefore solve this to any desired accuracy by 
an iterative procedure, in which the right-hand side is treated as a "source term" constructed 
using the solution for the B (x) obtained at the previous iteration. Thus we solve 

B (i) (x) - B^(x + 1) = a{x) B^ix + 2) , (3.7) 

where 

A 4 

a{x) = "(( 2; + l)2_ a 2)(( x + 2)2_ a 2 ) ( 3 - 8 ) 



and B m (x) = 1. The solution to ( |3.6[) is then given by 

B(x) = J2 B (l \x) = 1 + B w (x) + B {2) {x) + ■■■ . (3.9) 

i>0 

Each successive term in this sum is smaller than its predecessor by a multiplicative factor of 
A 4 , the small expansion parameter of the iterative solution{|] It should be emphasised that 
in the case we are considering here for the dyonic string, unlike the D3-brane discussed in 



[25], the constant a in (2/7) is itself A dependent. However, this does not alter the validity 
of the iterative procedure. 

It is straightforward to see from the above discussion that the solutions for the quantities 
B^ are given by 

B w (x) = 5^a(s+p) , 

p>0 

B m (x) = ^ a(x + px) a(x + pi + p 2 ) , 

pi>0p 2 >2 

B^{x) = Y Y a{x + Pl ) a{x + Pl + p 2 ) a(x + Pl + p 2 + p 3 ) , (3.10) 

Pl>0p2>2p 3 >2 



B (l) (x) = Y Y " ' Y a ( x a ( x + P 1 +-P2) • • • a{x + pi H Vpi) ■ 

pi>0p 2 >2 p»>2 

Defining 

oo q 

S(m ,mi,...,m q ) = ^ Yl {a(x + i + j)) mj , (3.11) 
i=0 j=0 

it can then be straightforwardly shown that the coefficients 5 w (x) are given by 

B w (x) = 5(1) , B {2 \x) = ±5(1) 2 - ±5(2) -5(1,1) , etc.. (3.12) 
1 There is a subtlety concerning this point, which we shall address in detail later. For now, we just remark 
that the statement is true provided that B (z> (x) is evaluated for a sufficiently large (positive) value of x, and 
that we can arrange for this criterion to be met. 
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The sums S(mo, mi, . . . , m q ) can all be performed explicitly, with the results expressed 
in terms of the Digamma and Polygamma functions. The expressions for 5(1), 5(2) and 
5(1, 1) were given in |25|| . We shall just present the expression for 5(1) here, since it is the 
only sum needed for our present purposes: 

qm _ (2x + 3)A 4 V»(l + a + x)-^(l-a + x) x4 

^ Uj "(4a 2 -l)(a 2 -(x + l) 2 ) + a(4a 2 _i) ' {6A6) 

where ip{x) = T'(x)/T(x) is the Digamma function. (Our definition of S(mo,mi, . . . ,m q ) 
includes a factor of (— A 4 ) n relative to the one in p5fl , where n = Y11=o m i ls the rank of 
S(niQ,mi, . . . ,m q ). Thus our 5(1) is (—A 4 ) times the one in p5| ; 5(2) and 5(1,1) have 
additional factors of A 8 , and so on.) As we shall discuss below, the series in ( |3.1lD defining 
5(1) is convergent, and free of singularities, provided that the denominators in a(x), given by 
( |3.S] ), do not happen to vanish for some particular terms in the summation. In particular, 
this means that for generic x the expression ( |3.13 ) should be finite for all values of the 



parameter a. There are apparent divergences at a = ±|, but these are actually spurious, 
and by taking the limit a — > ±| with due care, the expression (3.13) can be seen to be 
perfectly finite.^ 

3.2 The Floquet exponents 

The constant fj,, related to the Floquet exponent 2/i, can be determined from the observation 
that ^f(z) with exponent —2^ and ^>{— z) with exponent 2[i are two Floquet solutions that 



must be proportional to one another. From (3T) one then reads off that the coefficients 
C(n + fj) must satisfy the relation 

C[n + n) = pC(-n-n) , for all n , (3.14) 



where is an n- independent constant. In fact it suffices to impose ( 3.14 ) for just two 



consecutive values of n, say n = and n = 1, since all the rest can then be deduced using 



|3.2j ). Thus the content of ( 3.14 ) is just two independent equations, which determine the 



two unknowns (3 and \x. In particular, we may conclude that [i is determined by the equation 

C( M )C(-^ + l) _ 

c(m-i)c(-^) - • (3 - 15) 



2 A similar phenomenon arises also in the expressions that one obtains for the higher functions 
S(mo,mi, . . . ,m q ). For example the expression for the sum S(l, 1), given in pEj ]. has apparent diver- 
gences not only at a = ±4, but also at a = ±1. Again, these are spurious, and a careful limiting procedure 
reveals that S(l, 1) is indeed perfectly finite at these values. 
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The recursion relation ( fj^ ) can be re-expressed as G(x) = (V(x) — G(x + l)) -1 , where 
G{x) = C(x)/C(x — 1) and V(x) = (x 2 — q 2 )/A 2 . This may be solved to give the continued 
fraction 

G(x) = 1 x • (3-16) 

V[X) v(x + i)- V{x+ ^_... 



It follows that (|3.15[) , which can be rewritten as G(fi) G(l — yu) = 1, can be solved to give 



H as a power series in A: 

H = Y, A 2n = Mo + mi A 2 + 1X2 A 4 + • • • . (3.17) 

n>0 

Note that in our case we obtain a power series in A 2 , by contrast to the power series in 
A that arose for the D3-brane in [25]. This is because our coefficient a 2 in the Mathieu 



equation for the dyonic string is itself A dependent, as given in fl2.8| ). For the first four 
partial waves, t = 0, 1, 2, 3, we find, up to order A 8 : 

l-n. „_1 A /A2 7x2 A(2A 2 -3) X 4 (2A 2 -3)(8A 4 -12A 2 +5) ,6 
£_U. ^___ V A -1A ~ 2VA2 _ 1 A 8( A 2 -i)3/ 2 A 



A(240A 8 -1080A 6 +1798A 4 -1316A 2 +355) ,8 , /n/xKh 
48(A 2 -1)V2 A + U ( A )> 



1- ,, — 1 1a \2 3A 2 +8 \4 9A 4 +88A 2 ~36 \6 
1. /i-l-2^A 144A A 



135A 6 +3248A 4 -2032A 2 +576 \8 | 0( A 1O) (3 18) 



n. ,, _ 3 1a \2 4A 2 +9 \4 A(16A 2 +117) x 6 1600A 4 +25380A 2 +18387 \8 , <*y \ Hh 
l - M-2~3 ZiA T08~ A 1944 A 69984?) A + <A A ) 



— 9 . „ _ 9 _ 1 A \ 2 _ 15A 2 +32 X 4 _ A(225A 2 +1504) x 6 
t — O. [X — Z 4 ^AA 96Q A 115200 A 

16875A 4 +235712A 2 +140288 \8 , (n(\W\ 
55296000 V / • 

The general formula, up to order A 8 , turns out to be 

— 1(P-L~I\ A \2 2+1(1+2) (A 2 +2) a 4(l+l) 2 (3l 2 +6l+2)A+2l 2 (l+2) 2 A 3 x6 

M - 2 "I" t + i A l(£+lY (t+2) A I 1 (£+l) b (i+2) 2 A 

$ ( ,^*„^ {W + 1) 4 (15^ 4 + 60£ 3 + 55£ 2 - 101 - 12) 



' (e-i)e 3 {t+iy (e+2) 3 (e+3) ' 

+A{£ -1)(£+ lf{£ + 3)(15£ 4 + 60£ 3 + 80£ 2 + A0£ + 8) A 2 

+5(£ - l)£ 3 [l + 2) 3 [l + 3) A 4 } + 0(A 10 ) , (3.19) 



when £ = 2,3,4, .... In fact if £ takes a generic non-integer value, ( 3.19j ) is universally 



valid. The integer values for £ have to be treated carefully because the general solution for 
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ji develops poles at the integer values, and one has to resort to a case-by-case analysis. It 
seems that in fact only i = and I = 1 are not described by the general formula (|3.19| ).P| 



3.3 Asymptotic behaviour and absorption probabilities 

A remarkable property of the Floquet solution of the Mathieu equation is that one can 



re-express the series expansion (3.1) in terms of Bessel functions rather than exponentials, 
with the same coefficients C(n + jj) (up to (—1)" factors). (The explanation can be traced 
back to the fact that exponential functions and Bessel functions are closely related via the 



Laplace transform [32].) The required solution of the Mathieu equation, which is in-going 
on the horizon at z = oo, can be written in the form 

9(z)=MM= £ C c?f Jn(2Xe- z )H^(2\e z ) , (3.20) 

n=— oo " ' 

(The second solution is obtained by choosing an independent Bessel function, such as #4+2it> 
in place of the first Hankel function H^ 2 u-) 



The series (|3.9|) for ^f(z) is dominated by the n = term near the horizon, where z —* oo, 
and this allows the asymptotic behaviour in the z — > oo region to be read off easily. To 
determine the behaviour at spatial infinity, corresponding to z — > — oo, one can make use of 
the fact, mentioned earlier, that ^>{—z) with Floquet exponent 2\i is proportional to ^(z) 
with Floquet exponent — 2fi. In fact from ( |3.14| ) one obtains the result that 

M(-^z) = ^t^M^,-z) . (3.21) 



After some further steps described in detail in p5[ , in which one extracts the ratio between 
the amplitudes of the ingoing and outgoing waves at spatial infinity, one concludes that the 



3 In this section, we shall present explicit results for corrections up to and including 0(X 6 ) with respect to 
the leading-order results for absorption probabilities. One might think, therefore, that it would be sufficient 
to determine \i up to 0(A 6 ), since it has a A-independent constant term for its leading order. However, as 



will be seen in the calculation of the absorption probability P in ( 3.22 ) and (|3.23| ) below, there is a factor 
of \e 2l ^ 7r — e — 21 M7r | 2 in the expression for P, and so the fact that the constant terms in fi are integers or 
half-integers means that \e 2l ^ 7r — e _2l '" r | is a quantity of 0(A 2 ), rather than the naively-expected O(l). 
Consequently, in the calculation of \e 2l ^ 7r — e~ 2,f * n \ it is necessary to work with /i up to 0(\ a ), for the 
purposes of eventually obtaining corrections up to C(A 6 ) in the absorption probabilities. Similarly, the 
occurrence of a pole in r(— fi — a + 1 ) in (3.4) means that one must also work with fj, to order A 8 in 
that expression too, when calculating C(—[i) using (3.5). By the same token, in the D3-brane calculations 
described in |2q] , one should calculate fi up to 0(X lb ) in order to obtain corrections up to C(A 12 ) in the 
absorption probabilities. 
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absorption probability is given by 



P =] =WT? m> ( 3 - 22 ) 

\v~v L r + \x-x I 



where 



" = e2 ""- *=w- (3 ' 23) 

Note that when determining the orders in A to which one must work so as to obtain cor- 
rections to P to a specified order in A, one must be careful of a number of points. One 
such point was discussed in footnote 3, relating to the fact that the leading-order terms in 
ji are always integers or half- integers. Another point is that if the leading-order term in \ 
is a A-independent constant, then if this constant happens to equal ±1 then the quantity 
\x — X" 1 )] which would naively have been of 0(1), would actually instead be of higher 
order in A. This in turn would necessitate evaluating x to a higher order in A than one 
might naively have expected, for the purposes of finding the corrections to the absorption 
probability up to some particular order in A. This particular subtlety does not arise in our 
calculations in this section, but it will play a role in the calculations in section 4. 

The computation of the absorption probability is now reduced to the problem of calcu- 
lating the ratio C(-fx)/C(fi), where C(x) is determined from Q, Q, fl3~TI|) and fljmp . 
It is at this point that we meet a subtlety, which was not explicitly encountered in the gen- 



eral discussion in [26]. It arises because the parameters that occur in our Mathieu equation 
are not generic, but are instead rather specific ones that are fortuitously inconvenient. 

It is evident from ( |3.8| ) that the sums in ( [3.11 ) are rapidly convergent, and that the 
possibility of divergent behaviour would arise only if the denominators in a{x) were to go 
to zero for some particular low term in the summations. Since x, and a, are functions 
of the small parameter A, the issue that concerns us here is therefore not that the sums 
S(rriQ,mi, . . . ,m q ) might be infinite, but rather that they might acquire inverse powers of 
A that could counteract the manifest A 4 numerator in ( |3.8| ). If this were to occur, then 
our implicit assumption that each term B^ix) in ( |3.S| ) were smaller than its predecessor 
B( l ~^(x) by a factor of A 4 would be invalid. At the very least this would complicate 
the discussion because it would mean that we would have to work to higher order in the 
iterative solutions B M (x) than naively appeared necessary. At worst, it could be difficult 
to determine whether one had included sufficiently many terms in the iterative solution to 
give an accurate result to the desired order in A. 

First, let us note that no difficulty of this kind arises in the calculation of C(fi). This 
is because n = ±(£+ 1) + C(A 2 ), and a = \{l + 1) + C(A 2 ), and hence from it follows 
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that all the quantities a(x + i + j')/A 4 appearing in the sums fl3.11|) are strictly finite as A 
tends to zero. If we define the rank n = J2i=o m i f° r the sum S(mo,mi, . . . ,m q ), then it 
follows that the quantities A~ 4n S(mo, mi, . . . , m q ), evaluated at x = are all of the form 
of a non-zero constant +0(A 2 ) as A tends to zero, and so the iterative solutions i? w («) do 
indeed have the "naively expected" form i? w (u) = 0(\ ). 

The potential difficulties arise when one considers C(— /i), since now we see from ( |3.8|) 
that for sufficiently small values of the summation variables in ( 3.11| ) one will encounter 
denominators in a(x) that tend to zero as A tends to zero. This means that these sums 
S(niQ,mi, . . . ,m q ), evaluated at x = —fi, and the associated iterative solutions B (l) (—fi), 
will have A dependence whose leading order is at a more dominant power than the naively 
expected one. 

The easiest way to handle this problem is to make use once again of the recursion relation 
( |3.6| ) for B{x). By iterating it n times, we can express £>(— u) in terms of B{n — u) and 
B(n + 1 — fi). Thus we can choose re = £ + 1, so that the leading-order constant term in the 
A expansion of n — /i is precisely the same positive constant as occurs in the A expansion of 
re itself. By this means, we reduce the problem of calculating B(— u), and hence C(—fi), to 
the problem of calculating B(\(£+l) + 0(\ 2 )) and B(\(l+l) + l + 0(\ 2 )). By our previous 
discussion, these cases are easily evaluated, in terms of quantities 5(mo,rrei, . . . ,m q ) that 
have precisely the "expected" A 4fc leading-order dependences, where k = J2i m i- Note that, 
as mentioned in footnote 3, one must be careful to include one additional order in A 2 in 
the expansion for u, over and above the naively expected one, on account of the pole at 



A = in the quantity T(—fi — a + 1) that arises in the computation of C{— u) using (3.4) 
and Q. 

The upshot from the above discussion is that provided we handle the C(—fi) calculation 
carefully, it is necessary to include only the first iterative correction B w (x) = S(l) in order 
to obtain our results for absorption probabilities with corrections up to 0(A 6 ) to the leading- 
order result. (By the same token, the results up to C(A 12 ) corrections for the D3-brane, 
discussed in j25|, require the inclusion only of the first three iterations B^(x), B^(x) and 
-B (3) (x).) Of course the fact that each term B^(x) comes with an overall leading-order 
A 4i factor means that one need calculate the associated sums S(mo,mi, . . . ,m q ) only up 
to A n_4i corrections in order to obtain to obtain a results for B(x) to order re. (Assuming, 
as usual, that the necessary recursion using (|3.6|) is performed first, so that the B^(x) are 
evaluated at arguments x > a.) 

We find that the absorption probability for the ^'th partial wave has the general 
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structural 

47r 2 A 4+4£ n 

p -MFWTF|S l - l2 " (bgl) '' (3 ' 24) 

where A = e 7 A, and 7 is Euler's constant. The prefactor is chosen so that 60,0 = 1- Our 
results for the coefficients b n ^k with k < n < 3, for the first four partial waves I = 0, 1, 2, 3, 
are as follows: 

£ = : 6i )0 = 4A , 61,1 = -8A , 

6 2i0 = ±(32A 2 - 7) - f (A 2 + 2) vr 2 , 6 2 ,i = -4(10A 2 + 1) , 6 2 ,2 = 16(2A 2 + 1) , 

63.0 = A (64A 2 - 31) - f A (6A 2 + 17)vr 2 - f A(A 2 - 1) C(3) , 

6 3 .1 = -8A (22A 2 + 1) + f A (A 2 + 5) vr 2 , 

63.2 = 32A (6A 2 + 5) , 6 3 ,3 = -^fA(A 2 + 2) , 



t = 1 : 61,0 = 5A , bi,i = -4A , 

b 2 ,o = |(144A 2 + 14) - ±A 2 vr 2 , b 2 ,i = -±(63A 2 + 8) , 6 2 , 2 = 8A 2 , 
6 3i0 = ^A (2241A 2 + 704) - ^A(33A 2 + 8) vr 2 - |A 3 C(3) , 
6 3 ,1 = -iA (1251A 2 + 440) + § A 3 vr 2 , 

b 3t2 = |A (33A 2 + 8) , 63,3 = -f A 3 , (3.25) 



( = 2: 6 1)0 = fA, 61,1 = -§A , 

62,0 = gig (7472A 2 + 477) - ^A 2 vr 2 , b 2 ,i = -^(164A 2 + 9) , 6 2 , 2 = f A 2 , 

63.0 = 5^2 A (132992A 2 + 25443) - | f A(28A 2 + 3) tt 2 - § A 3 C(3) , 

63.1 = - ^3 A (3904A 2 + 657) + § A 3 tt 2 , 

63.2 = §A (28A 2 + 3) , 6 3i3 = -^A 3 , 



t = 3 : 61,0 = %A , = -2A , 

b 2 ,o = M 3875A2 + 171 ) " n A2?r2 » 6 2,i = -T20 (955 A 2 +32) , 6 2>2 = 2A 2 , 
b 3fi = 43^00 A (6053125A 2 + 786432) - T ^ fJ A(485A 2 + 32) vr 2 - ±A 3 C(3) , 
hi = "Tiioo A (255275A 2 + 27488) + \ A 3 tt 2 , 

^3,2 = ggA (485A 2 + 32) , 63,3 = -|A 3 , 

4 The absorption cross-section 01 is related to the absorption probability Pe by Ell: ae — + l) 2 Pi . 



13 



Note that our result for b\ t \ at I = agrees in magnitude, but not in sign, with the 
result obtained in 

For general values of £ it is quite straightforward to obtain expressions for the coefficients 
6i t o and &x,l m ( ft-24 ), using procedures analogous to those that we described previously. 



The only new feature is that in order to calculate £?(—//) in terms of B(x) with appropriately 
positive arguments (£ + 1 — fj) and (£ + 2 — fj,), one must now apply the recursion relation 
( |3.6| ) a total of {£ + 1) times. Up to order A 2 , it is easy to see that one need only solve 
( |3.6| ) to linear order in a(x), and evaluate B{£ + 1 — /i) and B(£ + 2 — ^) to zero'th order, 
B{1 + 1 - ii) = 1 + C(A 4 ) and S(l + 2 - //) = 1 + C(A 4 ), giving 

B(-fj.) = 1 + Y, a (-V + *) + ••■ • ( 3 - 26 ) 

i=0 

Furthermore, it is clear that only those terms for which the denominator in a(x) in ( |3.8|) 
has zeros at A = 0, namely i = £ and i = £ — 1, will give C(A 2 ) contributions. This leads to 
the result 

After further straightforward calculations, we find that the absorption probability is given 
by 

p '=(4^l 1 ^ l!kslt (W( 1+,(,+1, * tl) ) t '''l' M 

where i/j(x) = T'(x)/T(x) is the Digamma function. Note that in this case we have not 
absorbed the Euler gamma constant 7 in a rescaling of A in the logarithm. However, for 
all integer £ the Digamma function has the form ip{£ + 1) = R — 7, where R is a rational 
number, and so the usual A = e 7 A rescaling will enable the Euler constant to be absorbed 
in a natural fashion. 

4 Fixed-energy absorption, and "phase transitions" 

So far, we have calculated absorption probabilities in the regime where A is fixed, and A is 



a small parameter (see eqn. (2.8) for definitions of these quantities). These results describe 
the probabilities as continuous functions of low-energy incident waves, in a given dyonic 
string background. In this section, we shall study a different regime, where a, defined in 
(|2.8| ), is fixed. Namely, since the background depends on two charge parameters Q e and 
Q m , one can adjust the dyonic background by tuning A in such a way that a is fixed (for 
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a fixed value of A). This implies that 

\{l + l) 2 -a 2 = \ 2 A = \uj 2 M (4.1) 

is fixed, where M = Q e + Q m is the mass per unit length of the string. The low energy 
expansion with A <C 1 implies that A > 1 ( 1/A can be viewed as the expansion parameter). 
Note that in this case u> 2 M can be large, provided that A is sufficiently large. 

In order not to cause confusion with the previous section, we write the Mathieu equation 
as follows: 

+ (8A 2 cosh(2z) - 4a 2 ) * = , (4.2) 



where 



j9 i(^ + l) 2 -a 2 , s 

A 2 = ^ i- . (4.3) 

A 



In other words A is now viewed as a function of A, having first specified a. It is arranged 
so that A is small by choosing the charges so that A is sufficiently large. This should be 
contrasted with the viewpoint in section 3, where the parameter A is chosen to be small 
by taking u to be small, with the charges unconstrained and with a a derived quantity 
following from the expression in equation ( p.8| ). 

This equation is analogous to the Mathieu equation discussed in [^E| , in that a 2 is now a 
given fixed parameter, independent of the expansion parameter A. However, in our case a 2 



can be arbitrary, whereas in [25] the analogous parameter could take only discrete integer 
or half-integer values. We see from (|2.8|) that depending on the value of A, we can have a 2 
either positive, negative or zero. 

4.1 Generic a 2 > 

The Floquet exponents in this case turn out to be given by 
A 4 (60a 4 - 35a 2 + 2) A 8 



/U = a 



a(4a 2 - 1) 4a 3 (a 2 - l)(4a 2 - l) 3 
(6720a 10 - 18480a 8 + 15260a 6 - 4705a 4 + 413a 2 - 18) A 12 n( ^ 

4a 5 (4a 2 - l) 5 (a 2 - l) 2 (4a 2 - 9) ( ' ' ' 

This result will be valid for all generic values of a, but breaks down at the poles occurring 
when a is an integer or half-integer. We shall return to a discussion of these special cases 
later, and concentrate for now on the generic situation when ( |4.4j ) is valid. It is useful to 
define 

a = \L+l, (4.5) 
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so that for now, we shall consider cases where L is not an integer. For such generic values 
of the parameter a there are no subtleties involving the occurrence of poles or zeros in the 
calculation of C{— /i), nor in the expression \rj — rj~ 1 \ appearing in ( gjg ). We shall just 
present results for the first two corrections to the leading-order absorption. One has to 
distinguish between two cases, according to whether a > ^ or a < | (we shall consider the 
case a = | later). The reason for this distinction is that when a > ^ the next-to-leading 
order corrections are always of the form A 4 and A 4 log A, whilst if a < \ the next-to-leading 
order correction is of the form A 8a . For a > ^, we find that the absorption probability is 
given by 



Pe(a) 



7T 



: A 8a 



r(-2a) 2 -T{2af 



(l + 6i,iA 4 logA + 6 lj0 A 4 + - 



(4.6) 



where 



&i,o 



8 


r(-2a) 2 + T(2a) 2 " 




a (4a 2 


-1) 


r(-2a) 2 -T(2a) 2 



47r cot 27ra 



u , / 4a 2 



T -V(-2a)-^(2a)) , 



(4.7) 



a (4a 2 

where i/j is again the Digamma function. 

When a < ^, we find that up to the first correction to the leading order, the absorption 
probability is 



Pe(a) 



4A 8a T(-2a) 2 sin 2 2na ^ ^ T(-2a) 2 



cos Ana + ■ 



(4.8) 



r(2a) 2 1 r(2a) 

The next correction will be either A 16a , or A 4 log A, depending on whether a < \ or a > 7. 



4.2 



a 



|L + 1 



In the previous subsection, we saw that for generic values of a = + 1, corresponding to 
L non-integer, the leading-order absorption probability is of the form 



Pr. ~ A 



1L 



(4.9) 



This result is not applicable if L is an integer. This can be seen both from the resulting poles 
in the expression (|4.4|) for /i, and the singular behaviour of T(—2a) and the trigonometric 



functions in (4.6) and (4.8). In fact the absorption probabilities for the special cases L 



integer are precisely those calculated in [25|. These have the leading-order form 

P L ~ A 8+4i . 



(4.10) 
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To be precise, the absorption probabilities in these cases are given by ||2 

A 2 A8+4L n 

where A = e 7 A. The explicit results for the coefficients b n ^k with k < n < 3 were given for 
L = 0, 1, 2 in |2Ej. In order to calculate the corrections up to this A 12 order, it is necessary 
to calculate the Floquet exponent 1\i up to order A 16 , for the reasons that we discussed in 
footnote 3. They are given by 

r_n- ,, - 1 iy^ A 4 , 71 A 8 , 118511 a 12 12431571 a 16 , (nf \20\ 
L _U. ^ — 6" 216^5 62208^/5 " 55987200^ A + U { A ) ' 

r _ 1 . ,. _ 3 1 a4 , 133 a8 , 311 a 12 , 908339 a 16 , (n( \20\ (A^O\ 
L-l. [I - ^ - To A + 8640 + 3110400 A + 1254113280 A + U ( A )> \^ L/ ) 



2 12 11 ^ 8640 11 t 3H0400 11 ^ 1254113280 

r_o. _ o L A 4 _ 137 a8 , 305843 a 12 _ 64347197 a 16 , /n/ a20\ 

^ z • ^ z 30 54000 " r 1360800000 1714608000000 r ^ / 



It is also necessary to calculate the sums S(mo,mi, . . . ,m q ) up to rank n = J2i m i = 3. 
The expressions for n = 1 and n = 2 are given in [25]. 

There are in fact special values of the integer L which themselves must be treated as 
special cases. One of these is L = —2, corresponding to a 2 = 0. The other is L = — 1, 
corresponding to a 2 = These two cases are discussed in the next two subsections. 

4.3 a 2 = \ 

In this case, the Floquet exponent 2/i is given by 

fi= i +i A 2 -!fAA 6 + +^giA lo + 0(A 14 ) . (4.13) 

We find that the absorption probability is then given by 

P = 4vrA 4 (l + 16A 4 (logA) 2 -4A 4 log A - (| + § tt 2 ) A 4 + • • • ) . (4.14) 

In fact this result can also be read off directly from the t = expressions in ( 3.25Q , by 
setting A = 0. 



4.4 a 2 = 



When a = 0, we find that the Floquet exponent is given by 



» = V2A 2 (1 + f A 4 + -j-j^pA 8 + ^^ 2 + 0(A 12 )) . (4.15) 



We then obtain the absorption probability 



p = ^tct t 1 - f A4 <log A)2 - f (4C(3) - 3) vrmw + 0(A8) ) • «" 6) 
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where A = e 7 A. Note that as A — > the absorption probability ( 4.16 ) approaches zero not 
as a positive (integer) power of A, but instead as (log A) -2 . This result indicates that on 
the field theory side, the infrared limit with the specific choice of the dyonic background 
that ensures a = exhibits new features. For example, the two-point correlation function 
for the minimally-coupled scalar Q has a leading non-zero contribution that is inherently 
quantum in nature (as signalled by the log A dependence). Note that the a 2 = case 
that we have considered here is an example of the kind we mentioned in section 3.2, where 
the A- independent term in the expansion for \ ls equal to 1, implying that one loses the 
naively-expected leading-order constant when one evaluates \x — X~ 1 \- 

It would be interesting to explore the effects of infinitesimal deviations from the limit 
a = 0. It turns out that if one takes a special slice a 2 = 6 2 A 4 , where b = constant, one 
obtains a systematic expansion of \i around small A: 



/ f,2 25 , 9 l2 

M = A 2 ( Jl + b - + -^=A 4 + 0(A 8 )) . (4.17) 
V " V 1 + T 

In this case the absorption probability has the following form: 

P = H + (£«A). <^-l [C(3)(8^)-6] ,4p°^ A)2 +0(A«)) , 

(4.18) 

Note that as one approaches a 2 ~ 6 2 A 4 , only the sub- leading term in the A expansion 
depends on b 2 . 

4.5 a 2 < 

One can see from ( |4.3|) that it is possible to choose the parameters so that a 2 becomes 
negative, while keeping A small. In such a case, we may write 

a = i(3, (4.19) 

where /3 is real, and related to A by 

A2 _ & + ir + p _ 



In this case the Floquet exponents are given by (4.4) with a = i/3. The absorption prob- 
ability is now qualitatively different in form. In particular, the leading-order contribution 
becomes oscillatory as a function of A; we find that it is given by 

sinh 2 2vr/3 



sinh 2 2vr/3 + sin 2 (0 - 4/3 log A) 



(4.21) 
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where 

The oscillations become insignificant if >> 1. Again this region exhibits novel func- 
tional dependence of the absorpiton probability on A, thus signalling novel field theoretical 
phenomena on the dual side. 



5 Conclusions 

In this paper, we have studied the solutions of the minimally coupled scalar wave equation 
in the background of a six-dimensional extremal dyonic string. Owing to the fact that 
the equation is equivalent to the modified Mathieu equation, the wave-functions can be 
constructed exactly, in the sense that one can obtain exact solutions as power series in the 
parameter A 2 appearing in the Mathieu equation ( |2.7| ). This provides a systematic power- 
series expansion that is quite different in nature from a power-series solution of a differential 
equation in some restricted range of the independent variable. In particular, it allows one 
to have complete control over the relation between the small-distance and large-distance 
asymptotics of the wave functions; it is this information that is needed in order to calculate 
the absorption probabilities. 

Since the dyonic string background depends on two charge parameters Q e and Q m , 
there are two independent dimensionless quantities that can play the role of ordering pa- 
rameters in the iterative solution of the Mathieu equation. One of these is the quantity 
A = jlo 2 yJQe Qm: which can be arranged to be small by taking the frequency u of the wave 
sufficiently small. The other parameter is 1/A = \JQ e Qml {Qe + Qm) < 1/2, which provides 
a measure of the ratio between the electric and magnetic charges Q e and Q m . In section 
3 we calculated the absorption probability at low frequency, with A fixed, i.e. with a fixed 
dyonic string background. In this case the leading-order constant term in the expansion 
for a can take only discrete values, namely the integers and half-integers \{l + 1), where 
t characterises the angular dependence of the ^'th partial wave. The low-frequency expan- 
sion arises in even powers of A along with the characteristic dependence on powers of log A 
(A = e 7 A, where 7 is Euler's constant). These results are analogous in structure to those 
for other extremal black-hole or p-brane backgrounds. For example, the occurrence of log A 
in sub-leading terms in the energy expansion for the absorption probability is analogous 
to what is seen in other examples, such as the M-branes and D3-brane p0| , [2f|. By the 
AdS/CFT correspondence, the results of section 3 should yield information on correlation 
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functions in strongly coupled two-dimensional quantum field theory perturbed away from 
the infrared conformally invariant fixed point. 

In section 4 we calculated absorption probabilities in a different parameter regime, 
namely where a is held fixed, with 1/A being used as the small order-parameter of the 
iterative solution of the Mathieu equation. This can be viewed as an expansion in which 
an incident wave of fixed frequency (satisfying A <C 1) is absorbed in the dyonic string 
background with a range of charge ratios, in a limit where A is sufficiently large, implying 
that the charge ratio is very large, In this case the value of a is fixed, but can be chosen 
to have an arbitrary value. We first obtained results for generic positive values of a 2 . The 
expressions for Floquet exponents and absorption probabilities become singular when a 
takes integer or half- integer values, and these cases have to be treated separately; this was 
discussed in section 4.2. 

A case of particular interest which itself lies outside the general discussion of special 
cases in section 4.2 is when a = 0, discussed in section 4.4. In this situation, we find that 
the absorption probabilities depend inversely on powers of log A, which seems to indicate 
a correspondence to an intrinsically quantum regime of the two-dimensional field theory. 
The a = point can be viewed as a "critical point", since the absorption probabilities for 
a 2 > and for a 2 < (discussed in section 4.5) are qualitatively different. When a 2 > 0, 
the absorption probability vanishes in the limit A — * and approaches this limit with a 
positive (multiple of 4) power of A. On the other hand when a 2 < 0, the leading-order 
behaviour becomes an oscillatory function of log A. We also studied the region close to the 
critical point a = 0, by taking a 2 = b 2 A 4 , with b constant. 

Owing to the fact that the boundary conformal field theory for the extremal dyonic 
string is two-dimensional, this example provides an especially manageable arena for making 
concrete tests of the AdS/CFT conjecture. The fact that on the supergravity side one has 
full control over the study of the scattering problem provides a starting point for addressing 
the corresponding strongly coupled non-critical two-dimensional field theory. 
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